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cn ■ Abstract 

^ : 

■ The 0{Nfa^) electroweak radiative corrections to the Thomson scatter- 

Q,^ , ing matrix element are calculated for a general renormalization scheme with 

\ massless fermions. All integrals can be evaluated exactly in dimensional regu- 

Q_i' larization which in several cases yields new results that are remarkably simple 

^ • in form. A number of stringent internal consistency checks are performed. 

Q \ The Z-^ mixing complicates the calculation of 1-particle reducible diagrams 

' considerably at this order and a general treatment of this problem is given. 

Conditions satisfied by 0{Nfa^) counterterms are derived and may be applied 

^ \ to other calculations at this order. 



1 Introduction 



A century ago the electron was discovered by J. J. Thomson. Today the scattering 
process that bears his name serves as a means to define the most precisely measured 
physical constant, the electromagnetic coupling constant, a. Thomson scattering is 
the scattering of a fermion, most conveniently taken to be an electron, off a single 
photon of vanishingly small energy. It is only for such a photon that both energy and 
momentum for the process can be simultaneously conserved. The cross-section for 
Thomson scattering is 

where m/ is the fermion's mass and Q/ is its charge. It thus provides direct access to 
the strength of the electromagnetic coupling, a. Nowadays a is determined from the 
quantum Hall effect or from the electron's anomalous magnetic moment, but Thomson 
scattering remains the prototypical process by which a is defined in particle physics. 

A general renormalizable model contains a number of free parameters that must 
be fixed by experimental input in order for the model to become predictive. For 
highest precision the best measured quantities are used. In the case of the Standard 
Model of electroweak interactions there are three free parameters of the bosonic sector 
that are normally fixed using a, G^, the muon decay constant, and Mz the mass of 
the Z° boson. Of these a is by far the best known quantity experimentally. Its 
use in high-energy calculations introduces contributions from hadronic effects that 
can be mitigated to some extent using dispersion relations applied to experimental 
data Still an hadronic uncertainty remains but this may be eliminated using 

one extra piece of experimental data as input [Q. a then plays a key role in all 
calculations of precision electroweak physics. 

In recent years great strides have been made in the calculation of 2-loop and 
higher-order Feynman diagrams and their application to electroweak physics. In the 
processes considered to date, the main effort has been directed at the computation 
of the Feynman diagrams whereas the renormalization has normally been a very 
minor aspect. In the calculation of the (9(a^m^/M^) corrections to the p-parameter 
|^-|TD[ there is a single topology containing counterterms. In some calculations ||rT| 
the MS renormalization scheme is used and the counterterms are taken care of by 
discarding divergent pieces of diagrams. Such an approach is dangerous as it removes 
the check of the cancellation of divergences between diagrams and counterterms. The 
calculations are then performed in more than one gauge to test for consistency. The 
full corrections to the anomalous magnetic moment of the muon have been 

calculated |jl2|,[l^. This process first appears at 0[a) and so 0{o?^ counterterms are 



not encountered and much of the complexity associated with 2-loop renormalization 
in the Standard Model is avoided. 

As progress continues in 2-loop calculations, confrontation with the full complexity 
of 2-loop renormalization is inevitable. It is our aim here to calculate the 2-loop 
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corrections to Thomson scattering in a general renornialization scheme. We will 
limit ourselves to corrections that contain an internal fermion loop with the fermions 
assumed to be massless. These will be referred to as 0{Nfa'^) corrections where 
Nf is the number of fermions. Because Nf uniquely tags these corrections, they 
form a separately gauge-invariant set and can be expected to be dominant because 
Nf is quite large. Although the 0{Nfa'^) corrections arc a somewhat reduced set, 
the full complexity of the 2-loop calculation is manifest and all the intricacies and 
new features of the 2-loop renormalization are present. The calculation will allow 
us to obtain conditions on the counterterms which, once obtained, can be used in 
other calculations of this order. The resulting cancellation of divergences provides a 
powerful check that is an alternative to calculating in several gauges. 

In the present calculation, one of the greatest challenges is organizational. The 
various contributions must be arranged so as to avoid double counting and should be 
grouped in some logically consistent fashion. Decisions have to be made as to whether 
to carry wavefunction counterterms in individual diagrams, and thereby work with 
finite Greens functions, even when they can be shown to cancel in the final matrix 
element. The Z-7 mixing adds considerably to the complexity of the calculation 
at this order. The presentation chosen here is an attempt to achieve the goals of 
consistency and clarity. Contributions that are clearly related, such as the 2-loop 
Z-j mixing and photon vertex corrections are treated together as far as possible. 

In section 2 our notation is explained. Section 3 explains the renormalization of 
the Standard model and derives the relevant counterterms valid at 2-loops. Section 4 
reviews the calculation of Thomson scattering at 0{a). Section 5 discusses the role of 
wavefunction renormalization and Ward identities in the calculation. It is shown that 
the 1-loop Ward identities must be explicitly imposed. In section 6 the calculation 
of Thomson scattering at 0{Nfa'^) is described with separate subsections devoted to 
the various classes of contributions. Finally section 7 derives conditions that must be 
satisfied by the 2-loop counterterms. 

2 Notation and Conventions 

In calculating radiative corrections to 0(iVjQ!^), expressions will be obtained that 
contain the product of two 1-loop contributions and it will thus be necessary to 
distinguish between 1-loop fermionic and bosonic corrections. The order and type of 
a correction will be indicated, where needed, by a superscript in parentheses. Thus 
SZ^^f^ indicates the 1-loop fermionic part of the counterterm SZ. The 1-loop bosonic 
corrections are denoted by the superscript ^^^"^ and the superscript indicates both 
together. The superscript ^'^^ when used here means the full 0{Nfa'^) correction. 

Ultraviolet (UV) divergences will be regulated by dimensional regularization in 
which n denotes the complex number of space-time dimensions. Most loop integrals 
will be given in exact form rather than expanding about n — A. In fact, keeping 
the full n dependence helps display some dramatic cancellations that occur between 
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different classes of Feynman diagrams. It is assumed that nowadays, with the wide 
availability of computer algebraic manipulation programs, the exact results are easily 
transformed into series expansions when required. With this in mind some expressions 
are conveniently and compactly written in terms ofe = 2 — n/2. 

Two-point functions for the vector bosons, n^,^(g^), can always be divided into 
transverse and longitudinal pieces. 

Here we will be exclusively concerned with the transverse part Uriq^)- The subscript 
V' will therefore be dropped. 

Throughout this work the Euclidean metric is used with the square of time-like 
being negative. The calculation is performed in 't Hooft-Feynman, R^=i, gauge. 

A fully anti-commuting Dirac 75 will be assumed. This could only lead to dif- 
ficulties in fermion loops that generate the antisymmetric e tensor, such as internal 
fermion triangles. In that case, however, when one sums over a complete generation, 
anomaly cancellation guarantees that additional terms cannot appear. 



3 Renormalization of the Standard Model 

3.1 Renormalization of the Bosonic Sector 

The bare lagrangian, £° is the true lagrangian of the theory. The renormalized 
lagrangian, and countcrtcrm lagrangian, 6C, satisfy £° = + SC. The bare 
lagrangian of the Standard Model for free neutral gauge bosons is 

^'w - - li9^K)9^Bl - {g'^Bl - g^W,,Y (2) 

where and B'^ are the bare SU{2)l isospin and U{1) hypercharge fields and 
and g'^ are the SU{2)i and f/(l) coupling constants respectively. (M^)'^ is the bare 
mass squared of the W boson and (M|)° will be used to denote that of the Z° boson. 
They are constructed from parameters of the Higgs sector from which the relation 

(M|)0 ^02 + g,02 ^ ) 

may be derived. The renormalized and countcrterm lagrangians are obtained by 
writing the bare fields, coupling constants and masses in terms of the corresponding 
renormalized quantities and their counterterms, 

^{\ + 8Z^)\w g'^g^h {Mlf = Ml^ml, (4) 

= + g'^'^g' + H {Mlf^Ml + ml (5) 
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These conventions differ from those used Ross and Taylor |Tj] and are more closely 
consistent with those of Aoki et al 

The weak mixing angle, Ow, is defined so as to diagonalize the mass matrix of the 
renormalized fields and hence the renormalized and B fields are then related to 
renormalized Z field and photon field, A, by 

W^= ceZ + seA, (6) 
B = -sgZ + ceA, (7) 

where sg and cq are sin 6w and cos Ow respectively. 

The transverse parts of the 2-point counterterms for the photon, Z^ and Z-7 
mixing are obtained by substituting eq.s(§) and into (||). Keeping only terms that 
can contribute up second order at = gives 

^aIa/XaaL = -]^Ml{2seCe6- + Uzz.f (8) 



z z 

.AAA/^X^AAA/^ = -SMl - M\bZz 



Ssece ^ V J 2 

^9 1 

+Sece6-{5ZB + 6Zw) - ^seced- — \ (9) 

9 J 



-445! - 2^^5.5 Zz^ + 2-5_(5Zb + 5Zw) - AsldJ-^ 
ce ce g 



where here and in what follows we define 



The relation 



(10) 



6Z^ = Cq6Zb + sl6Zw, 6Zz = s16Zb + cl6Zw, 5Zz^ = seCe{5Zw - 5Zb), (U) 
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M| ^^^^ 
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follows from eq.(|^) and is valid to first order. Note that at this order the photon 
develops a mass counterterm which is a reflection of the fact that the renormalized 
field, Afj^, is not the same as the physical photon. In principle one could rediagonalize 
the neutral boson mass matrix at each order by redefining 6w but it is cumbersome 
and unnecessary. 

In the corrections to Thomson scattering the charged VT-boson appears as an 
internal particle. Its 0{a) 2-point counterterm is 

JL^x(yJ^ = -6Z^^\q^ + M^)5^, + 5Z^^\^q, - SM^^S^, (14) 



Note that the first term in eq.(0) contains an inverse propagator, (g^ + M^r), 
and so for internal W^s it will cancel with wavefunction counterterms at vertices. Of 
course coupling constant counterterms and external particle wavefunction countert- 
erms still need to be included but internal W^s therefore effectively only generate 
2-point counterterms of the form dz'^q^q^ — 5M^^-'(5^,^. In practice this provides 
a very convenient way of eliminating much of the labour in calculating diagrams 
constructed by inserting Oi^NfO.) counterterms in 0{a) diagrams. 

In the on-shell renormalization scheme the 1-loop counterterms are 

6Mf^ = ReU^Xi-M^), (15) 
6Mf^ = Ren«(-M|), (16) 



se 


ngi(o) 




ce 




Se 


ngi(o) 




Ce 


Ml 



Sg'^'^ ^r^mur.. seii'z'io) i^/n«(-M2) ng)(-M|; 



g' 2 ^""^ ^ ce Ml 2 Mf 



z 



In the MS renormalization scheme the counterterms are just the divergent parts of 
these plus certain other constants. By direct calculation of the diagrams concerned 
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the divergent parts of the 1-loop counterterms are found to be 



r2(l&) _ (j^\M^ 
167r2 7 6c^ 



- (Tlr^)S(31-^-25)A (19) 



= -f#.l^(424-74.^ + 25)A (21) 



5^(16) / ^2 X 43 



(25) 



^4^^ = -2^^ = - ( ^ ) (26) 



g Vl67r2 y 12 



A (27) 



in which A = tt ^r(e). In all cases the fermionic counterterms have been summed 
over a single complete generation. 

3.2 Renormalization of the Fermionic Sector 
3.2.1 Fermion 2-point counterterm 

The bare lagrangian for a free fermion, /, is given in terms of the bare fermion field, 
■0°, and fermion mass, m° by 

£v, = -V'°(^ + m5)V^°. (28) 

In the Standard Model the left- and right-hehcity, '4>l and components of the bare 
field are renormalized independently. The rcnormalized fields and fermion mass are 
defined from the corresponding bare quantities by the rescalings 

iIjI = {1 + 5Zl)Hl. (l + 5Z^)^ij, mj = m/ + 5m/. (29) 
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W Z Y 



(a) (b) (c) 

Figure 1: Fermion self-energy diagrams 



This leads to a 1-loop counterterm 

— > = -i^{5ZLiL + ^Zbtir) - (^ ^^^ 2 '^^^ ) ~ ^'^^ ^^^"^ 

= -]^{5Zr'^l + 5ZLiR){ii) + mf) 

~]^{i^ + mf){5ZL-iL + ^Zr-hi) - 5mf (31) 

where the form (^Tj) is expressed in terms of inverse propagators and is useful for 
demonstrating the cancellation of fermion wavefunction counterterms where they oc- 
cur in internal loops. 

The finite pieces of the counterterms will depend on the particular renormalization 
scheme that has been chosen but the divergent part is common to all schemes. This 
can be found by computing the 1-loop Feynman diagrams contributing to the fermion 
self-energy. These diagrams are shown in Fig.|l| and give 

^,„,(^)Qj,.,„,),-._^r(2-|) (32) 

where terms that are suppressed by factors m^/M^ relative to the leading ones have 
been dropped and and [3Rf are the left- and right-handed couplings of the Z^ to 
the fermion, 

hf = '-^l^^, (3nf = (33) 
Ce Ce 
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The first term of eg. (|3^) comes from Fig.|l]a, tfie second and third terms from Fig.|l]b 
and the last two from the pure QED diagram, Fig.[^c. 

In the on-shell renormalization scheme, defined by setting the renormalized fermion 
mass rrif equal to the pole mass, the fermion mass counterterm is then given by 



ruf 



' ^ «}^i(-,)*-|:^r (2 - ^) (34) 



167r2 y (n-3) 



3.3 Fermion-Boson interaction lagrangian 

The bare interaction lagrangian between the neutral gauge bosons and fermions is 

= ^9%^hMwl + ^9"^iP,l,^lBl + ^9"^^h.KBl. (35) 

ip represents the fermion wavefunction. Here 7^ are the usual Dirac 7-matrices. The 
electric charge, Q, of given fermion flavour and helicity is related to its hypercharge, 
Y,hYQ = ts + YL = Yn. 

Substituting eq.s(§) and @ into (|35|) yields the vertex counterterms for the pho- 
ton. To second order in the counterterms for the coupling constants and bosonic fields 
and first order in the fermionic counterterms this is 

--5Z^ [^-5Z^ - -J ^5Zs [j5Zs - — 

+^gQsel^.lR\5ZR + f % + Uz,^ - Uzb (]sZb -%]] (36) 




9' 2 2 ^ V4 9' 

and for the Z^, although it is not required here, the vertex counterterm is 
Z 




,5Z, + t,ce f ^ + Uz^) - (g - Hfl (^4 + \^Zb) 
Vf2 J C0 \g' 2 J 



—^ce6Zw ( -SZw ^ I + ^^—-—^—SZb ( -SZb 



2"- - V4 9 J 2 ce ^ V4 ^ 9' 

(37) 



+^97^^7Ri /3rSZr - Q— — + -6Zb + —SZb -6Zb - 

{ Ce \g' 2 J 2 cg \4: g' 
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Here 7^ and '~)r are the left- and right-helicity projection operators. It is a simple mat- 
ter, using eq.(|3lD and (p6D, to show that the fermion wavefunction counterterms, 5Zl 
and 5Zr, cancel between vertex and 2-point counterterms in all Feynman diagrams 
of interest here. We will therefore not consider them further. 



3.4 Higgs field lagrangian 

The bare Higgs field after spontaneous symmetry breaking generates charged and 
neutral Goldstone scalars (0^)°, Defining the renormalized scalars fields and 
wavefunction counterterms via the relation 



+ (38) 



leads to 0{a) scalar and vector-scalar mixing counterterms 



= -5Z^q^-5l3 (39) 



.AA/^>^ 



w 



5Za, + -5Zw + — 
2 g 



(40) 



here 5/3 is a quadratically divergent 1-point counterterm. It will be assumed that it is 
adjusted to exactly cancel the tadpole contributions and can therefore be ignored. A 
detailed complete renormalization of the Higgs sector along the lines of ref. 
be found in ref. where it is shown that 
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can 



5Z. 



(1) 



5M, 



2(1) 

w 
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(41) 



Note that the fermionic part of this counterterm is finite and in the MS renormaliza- 
tion scheme is therefore set to zero. This is also true of the mixing between neutral 
scalars and vector bosons. 



3.5 Gauge-fixing lagrangian 

In the calculations performed here 't Hooft-Feynman, -Rg=i, gauge will be employed 
for which the gauge-fixing lagrangian is 

C,,. = {2 \d,W^ - iMw^^\^ + {d,Z^ - ^Mz<Pzf + {d.A^f] (42) 

with ^ = 1. At tree level this has the advantage that the mixing between vector 
bosons and Goldstone scalars is canceled between the vector-scalar interaction la- 
grangian and the gauge-fixing lagrangian, C-^s.- In order to satisfy Ward identities 
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£g.f. is constructed from renormalized fields. Mixing counterterms then reappear in 
2-loop diagrams thereby nulhfying an important advantage of i?^ gauges. Certain 



authors ||18| have chosen to replace the renormalized fields and masses in eq.(^) 
with bare one and satisfy the Ward identities by renormalizing the gauge parameter, 
^. Mixing counterterms are thus eliminated but no reduction in labour is achieved 
because gauge-parameter counterterms now appear and the two approaches are for- 
mally equivalent. The latter, however, goes against the notion of renormalization 
as a rescaling of the physical parameters. In the present work we follow Ross and 



Taylor [jl4| leaving £g.f. unrenormalized. 

As is well-known it is only correct to include £g.f. of eq.(^) if the correspond- 
ing Faddeev-Popov ghost lagrangian is included as well. It can be shown that the 
fermionic part of the 0{Nfa) 2-point ghost counterterms can only take the form 

r = -5Z^{q^ + M'^) (43) 



where is the ghost wavefunction counterterm. Note that there is no ghost mass 
counterterm per se which is consistent with the requirement that the ghosts appear 
only in closed loops and do not couple directly to fermions. The presence of the inverse 
propagator (g^ + M^) in the ghost counterterm means that contributions from this 
2-point counterterm cancel against diagrams containing ghost vertex counterterms 
and thus at O^NfO^) the ghosts effectively go uncorrected. 



4 Charge Renormalization at 0{a) 



The Thomson scattering amplitude has been calculated by a number of authors 19- 



|22| . The results of ref. p2[ are given for a general renormalization scheme assuming 




only that renormalization of the bare parameters of the model takes the form given 
in eq.(§) and eq.(^) and the same conventions are adopted here. For external photon 
momentum = 0, the sum of 1-loop photon- fermion vertex diagrams and external 
fermion leg corrections is 

= 4''^(0)7m7l = 2|^2t37.7L(vrM^)-T(e) (44) 

The corresponding vertex and external leg corrections for the Z^'-boson are 

= 4''^(0)7m7l = z^2t37.7L(vrM^)-T(e) (45) 

In both cases care has been taken to use the Feynman rules obtained from the renor- 
malized lagrangian without applying the relation, gse = q'cq. This is important when 
one comes to derive the 1-loop counterterm insertions at 0{Nfa'^). 
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The Z-7 mixing at = that contributes to the Thomson scattering matrix 
element is given by 

-^K>^ =5^.ng^)(0) = -2 + f ^ seCe6,.Ml{^Ml)-^T{t) (46) 

= -2-|^seCe(5^.M|(7rM^)-T(e) (47) 

where again in eq.(^) care has been take not to apply the relation gsg = g'cg. The 
1-loop fermionic corrections to the Z-7 mixing, ng^^(O), vanish at = 0. 

The photon self-energy is guaranteed to be purely transverse by gauge invariance 
and will be written 

-K>^ = iq%u - q,qu)%,iq') (48) 



with 



^ ':(-M^)-T(6) - '-^^{^M-^reT{e) (49) 



167r2 '12' 167r2'3' 
9'sl\ , 2 



167r2/ V 3 



3 + -e (7rM4.)-T(e) (50) 



and 



n;W'(0) = -fj^|'r-T(.)j:QK"'?)-' (51) 

/ 

where the sum is over internal fermions. When these internal fermions are quarks 
eq. (|5TD cannot be evaluated reliably because of strong QCD corrections. In that case 
one writes 

n;(/) (0) = Re n;(/) {f) - [Re n;(/) {<f) - n;(/) (o)] (52) 

with (f being chosen to be sufficiently large that perturbative QCD can be used. For 

Ig^l ^ m'j 

„V r(2-5)r(s")' 
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where we have summed over a complete fermion generation. The quantity in eg. 
in square brackets is obtained from the experimentally measured cross-section, (Jh{q^) 
for e+e~ — ^hadrons by means of the dispersion relation 



Ren;(/)(g^) - n;(/)(o) = / ^ /'V^. dq\ m 

47r^a — + le 

where is the mass of the vr*^. This dispersion integral has been evaluated most 
precisely for (f = — M| 

The 1-loop counterterm contributions to Thomson scattering may be obtained 
from eq.s(^) and (pG]). Combining all the contributions gives the result that to 0{a) 

VA^ = e^l + -n n(0) - + s,^ + c,^ j (55) 

in a general renormalization scheme. The quantity a that appears on the left-hand 



side of eq. (|55D is the experimentally measured value = 137.036..., and all pa- 
rameters on the right-hand side are the renormalized parameters in the particular 
renormalization scheme that has been chosen. 

Note that all dependence on the wavefunction renormalization counterterms, SZw 
and 6Zb has canceled and, as a consequence, one can safely set = 6Z^j^^ = 

as was done in ref.s |]TU|,|22|. This choice leads to divergent Green functions that, 
however, combine to yield finite expressions for physical quantities such as eq. (|55|) . 
This cancellation of divergences is a useful and stringent check. It will be seen, 
however, that at 0{Nfa^) the 0{Nfa) wavefunction counterterms must be explicitly 
included in order to obtain physically correct results. This amounts to imposing the 
1-loop Ward identities by force. 

5 Wave function Renormalization and Ward Iden- 
tities 

As noted in the foregoing section at 1-loop order, one has the option of setting 

= 54'^ = (56) 

because physical results such as eq. (^5)) are independent of them. The same cancel- 
lation of the dependence on 5Z'^ and 54^"* ^^at occurred at 1-loop will obviously 

(2) (2) 

occur for 5Zy/ and 5Z^ at 2- loops but then, as will be demonstrated, the condition 
(^) cannot be maintained and the 1-loop Ward identities must be imposed explicitly. 
Actually this feature is seen to be quite general. At (9(a") the wavefunction coun- 
terterms 5Z^'^'^ will cancel out in physical expressions but lower order wavefunction 
counterterms must be included. 
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(a) (b) (c) 







(d) (e) 

Figure 2: Diagrams containing a 1-loop fermionic counterterm and a factor 
(7rM2.)-T(6). 

Consider the diagrams shown in Fig.|^. The counterterms denoted by 'x' are the 
fermionic pieces only and the self-energy and vertex blobs represent bosonic radiative 
corrections. The blobs containing an 'x' denote the bosonic 1-loop diagrams with 
one-loop fermionic counterterm insertions. These diagrams are all proportional to 
the fermionic part of a 1-loop counterterm and the quantity (7rM^)~'^r(e) and there 
are no other diagrams having this functional dependence. The matrix element for 
Thomson scattering is expected to be proportional to the charge, Q, of the external 
fermion and independent of its weak isospin, t^. It follows that those parts of the 
diagrams in Fig.^] proportional to ^3 must cancel amongst themselves. By explicit 
calculation this is found to be 

(57) 

The 1-loop Ward identities require that 

+ ^ = (58) 

bzliP + 'j^ = (59) 
^ 9 

where eq.(p8D is true for both fermionic and bosonic counterterms separately. Hence 
(^) correctly vanishes provided SZy/ and SZ^ are included in a manner consistent 



with the 1-loop Ward identities. The conditions (^) and (|59D will therefore be applied 
where needed in the following. 

It will also be useful to note that in any renomalization scheme 

2 g \lo7r^ 



-5Z^'') + ^ + 2(-|-^)(7rM^)-T(6). = finite (60) 



2—— = finite (61) 
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The former vanishes in the on-shell renormahzation scheme and the latter in MS. 

Although rather arduous, we have checked that the sum of terms proportional 
to ts from pure counterterm contributions also vanishes. This is a useful check of 
combinatorics and the counterterms (|^) and (|36|) . In this case however the cancellation 
happens without having to impose the 1-loop Ward identities explicitly and is vahd 
beyond C»(iV/a2) and up to cfa^). 

It can also be shown that the 0{Nfa'^) corrections from one particle reducible 
(IPR) proportional to cancel amongst themselves. 



6 Charge Renormahzation at 0[Nfa^) 

All results given in this section will assume one massless fermion generation and in 
loops in Feynman diagrams will be summed over all fermions in that generation. 



6.1 One-particle reducible diagrams 



The presence of mixing between the and the photon greatly complicates the cal- 
culation particularly in the counterterm and one-particle reducible (IPR) sectors. 
Baulieu and Coquereaux [|l6l have shown how to treat Z-7 mixing to arbitrary order 
in a. Their results can be applied straightforwardly to neutral current processes such 
as e^C fi~^fi~ but it is not immediately clear how to treat the case of Thomson 
scattering where the photon is external. In ref. it was shown how to rearrange the 
expressions obtained by Baulieu and Coquereaux in a form that displays the exact 
factorization of the residue at the pole of a resonant matrix element that is known 
from S'-matrix theory to occur even in the presence of mixing. The same procedure 
can be used to obtain an exact expression for the residue at = for some neutral 
current process such as e'^e" — > . In that case the initial state residue factor is 
found to be 



nz-y(o) 



- UzziO) 



(62) 



which is precisely the Thomson scattering matrix element up to crossings. Here 
Viziq'^) and Vij{q'^) are the exact, all order, e~^e~Z and 6+6^7 vertex corrections 



respectively. The self-energy and mixing corrections Ilzz{(f), ^Z'yiQ^) and 11^7(0) 
are also the exact expressions. Expanding the square root generates the appropriate 
factors for IPR diagrams and the correctness of the procedure is confirmed by the 
cancellation between the IPR counterterm contributions proportional to with those 
coming from higher-order counterterms as described in section ^ 
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Figure 3: 0{Nfa^) one-particle reducible (IPR) diagrams. 

It was shown in the previous section that the wave function counterterms, SZ^^ 
and 6Z^j^^ must be included in a manner consistent with the Ward identities and as 
stated above it can also be shown that the contributions from the IPR diagrams pro- 
portional to ts cancel amongst themselves. The calculation can therefore be organized 
in such a way that the IPR diagrams and their associated counterterms together form 
a class that is separately finite and proportional only to the charge, Q, of the external 
fermion. This also means that there will be a separate cancellation of the divergences 
of (9(A^/a^) one-particle irreducible (IPI) diagrams with their associated countert- 
erms. In particular, it follows that the divergence structure of the IPI diagrams is 
not influenced by the IPR sector. 

The IPR self-energy diagrams contributing to the Thomson scattering matrix 
element are shown in Fig.^. The self-energy blobs are indicated to be fermionic or 
bosonic contributions by the 'f or 'b' below them. The associated combinatoric 
factors are also given. Note diagrams containing the 1-loop vertex corrections, ( ^4|) 
and (|i5D, do not appear because they are proportional to and have been shown to 
cancel as discussed in section |^. Discarding all terms proportional to ts from the IPR 
diagrams their contribution is 

+ ^94p'z'\0).Ilf/\0) + ^,.,^n;(f (0).n?f (0) (63) 



where 



n;«(o) = n;«(o) + 2(^. 



r Cn ; — 
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Figure 4: 0{Nfa'^) photon vertex corrections. 



The eq.s(0) and (|5^) have been used along with the fact H^^ (0) = 0. The 



quantity 11^^ (g^) is the derivative of 11^^ (g^) with respect to g^. It will have a 
hadronic component for = that may be obtained using methods described in ref. 



2J] . This hadronic contribution is distinct from the one that appears in n?^^^(0) that 



was discussed in section ^. For the leptons may be evaluated perturbatively 

from 



r-vr 

167r2 3 



Qfi 



m. 



(64) 



In the on-shell renormalization scheme all terms vanish identically because of the 
definitions of the 1-loop counterterms. This does not eliminate hadronic contribu- 
tions, however, because they will reappear when the counterterms obtained from 
charge renormalization are used in other calculations and will give rise to an hadronic 
uncertainty. 



6.2 Vertex and Z-'j corrections 
6.2.1 Diagrams 

Representative topologies for the OlNfa"^) photon vertex diagrams contributing to 
Thomson scattering are shown in Fig.^. These may be calculated using methods 



described in ref. |£5[. Diagrams of the type Fig.^-f containing virtual photons or 
Z^'s, instead of W bosons, cancel by Ward identities and we have explicitly checked 
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that this occurs. The sum of all diagrams in Fig.^ is 



16^ 8^3 W.^^ T{4-n)T[2--)T[^- 



(65) 



exactly for all n. 

Representative diagrams contributing to Z-7 mixing in 0{Nfa'^) are shown in 
Fig.^. Again methods for calculating the individual diagrams may be found in ref. 
Upon summing all diagrams together one obtains the remarkably simple result 

-^K>^ = 5/..ng)(o) 



(_|L)%.„.„,,„M|(«r,4-„,r(.-|)r(| 



(66) 

The 0{Nfa^) diagrams when added together must form a pure vector current 
proportional to the charge, Q, of the external fermion and independent of its weak 
isospin, ^3. Contributions from the photon self-energy, which will be dealt with later, 
are obviously of this form and cancellation of terms proportional to ^3 is expected 
between the vertex corrections, given above, and the Z-7 mixing when it is coupled 
to the external fermion. This is indeed borne out and one obtains 

9 . 2..ni-)(o) 



v|:^\Q)l,lL + ^-l,{h^L-slQ) 



2 \ 2 



6.2.2 Counterterm Insertions 

The 0{Nfa'^) corrections coming from the insertion of 1-loop counterterms into 1- 
loop diagrams may be calculated using the expressions for the 2-point counterterms 
given in section 3 and simplified using the Ward identity (|5D|). Note once again that 
contributions coming from the first term in the 2-point counterterm for the W boson, 
(0), cancel against vertex counterterms. The result is 

V^'\0h,jL = ^-^19 (yI^) 3t3Se7,7L(vrM^)-T(e) 

6M^'' . ( 9' 



W (yI^) t3SWi(7rM^)-^er(e) (6^ 
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Similarly the 0{Nfa'^) contribution from 1-loop counterterm insertions into the 
1-loop Z-7 mixing is 



^Mr^(|^)^(-M^)-(2-e)r(e) (69) 



Individual self-energy diagrams contributing to eq.(J6) contain pieces proportional 
to g'^SQCeMy^ and g^i^sl/ ce)M'y^ and there is a very complex and intricate interplay 
between diagrams containing W boson 2-point counterterms (p!^), scalar countert- 
erms ( p9D the vector-scalar mixing counterterms (^) to produce an overall result 
proportional to (7^(s6)/ce)M^. This, when connected to the external fermion by a 
propagator leads to a result proportional simply to g^se of the same form as ^^^^^''(0) 
in eq.(|68|). 

Combining eg . (1681) and eg . (1691) yields the total Oi^NfO^) contribution counterterm 
insertions in the vertex and Z-7 mixing, 



^i-y (0)7m7l + i—l^,[hlL - SqQ) ' 

Cg IVl ^ 

2(1/) \ / „2 



+ VTe^J Q.^7.(vrM^)-l2 - e)r(e) (70) 

Note that the finite terms, proportional to t3er(e) have canceled. The remaining 
part proportional to ^3 contributes to eq.(|^ and therefore cancels with other terms 
as discussed in section ^ It will be discarded and only the second and third terms in 
eg. (|70D , proportional to Q, will be retained. 

6.2.3 Counterterms 

The expressions for the counterterms given in eq.(H), (H) and (^) are correct to 0(«^) 
and need to be specialized to Oi^NfO^). Using the Ward identities, ( |55D and (^), the 
Z-7 mixing counterterm of eq.(y) becomes 



9 9 9' 9' 



"nil~U^^ ' 9 J 



+ ^ (71) 
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and the photon vertex counterterm, eg . (|36|) , yields 

'(1/) 




I V ^ 9 9 9 



-i9SeQlt,^ —■ — (72) 

9 9 

Their contribution to the Thomson scattering matrix element together is 
g ' ^ ^ ^ ~^ ^ ^9^hsei^iL 




^ f Jg^'"^ 6g^'f^ .5g'^^') 5g'^^f^\ 
-igseQ .?>\se—^.—^ + Cg—^.—^\-1^ (73) 

As discussed in section ^ the part proportional to can be shown to cancel against 
products of 1-loop counterterms coming from IPR diagrams and will therefore be 
discarded. 

It is convenient at this point to define a quantity, 11^^(0), obtained by combining 
the parts proportional to Q of eg. (|6^) , eg.(|70D and all but the last term in eg. ([73|) . 
Hence 

flfi(o) - (i^j8..c.M|<^r(4-„)r(2-|)rQ) 

+ ^ .,c,M|(7rM^)-T(6) 



g \l67r" 

' ^ i ^ ^ seCeMl{7,M^)-%2-e)V{e). (74) 



V 167r2 



6.3 The Photon Self-Energy 
6.3.1 Diagrams 

Representative topologies for the diagrams contributing to the photon self-energy, 
n!^^(0) of eg.(^) at Oi^NfO^) are shown in Fig.^. Calculation of the photon self- 
energy involves projecting out the transverse parts of individual diagrams followed 
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by differentiation with respect to the external momentum squared using techniques 



described in ref. |^5[. Individual diagrams are not separately transverse but we have 
checked that the longitudinal part vanishes when all diagrams are added together. 
The result is 

,2. \ 2 



x(!i^(.M|)-r(5-„)r(2-f)r(f-i 

f \ V / Tl /\/f2^--2■p/Q \r'/o 



f 



ghe \ ' 16(n - 2)^„_,^^^,^„_,^ _ r f2 - ^ 



levrV 3n ' V 2/ V 2 

,2 ( (^Lf + (^Rf 



+ (^V-^vr"-(M|)t-rf3-^)rf2-^ 

^ 167r2 y 3(n-2) ^ V 2/ V 2 

/ 

^ g^se \ ' 43^ (5^2 - 33n + 34) r /"o """iX^n 

+ ^16^ J T n(n - 5) " ^ (^^ - - j T (2 - - j )^ Q,,,., 



4 / 2 



(75) 

Contributions that are suppressed by factors m'j/M^ relative to the leading terms 
have been dropped. The first term on the right hand side of eq.(^) comes from 
diagrams Fig.^-e that contain an internal W boson, the second comes from diagrams 
Fig.^&g containing an internal Z^. The last term in eq.([75|) comes from diagrams 
Fig.pi&i that are pure QED in nature. Contributions for which the fermion mass 
can be safely set to zero without affecting the final result were be obtained using the 
methods described in ref. |^5||. The terms in which the fermion mass appears are 
obtained using the asymptotic expansion of ref. ||2^. It should be noted that setting 
771/ = in Fig.^-g does not immediately cause obvious problems in the computation 
because the diagram still contains one non-vanishing scale. A certain amount of care 
is thus required to identify situations in which the fermion mass cannot be discarded. 
In the case of the photon self-energy, the need to include such terms is indicated 
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by divergences proportional to the In mj in the counterterm insertion diagrams. All 
diagrams that yield m^-dependent terms can be split in two by a cut through two 
internal fermion propagators. The contributions are therefore precisely those that are 
accessed via dispersion relations (0). 

As it stands the eq.([75|) for 11^^^ ■'(0) contains a divergence, in its fourth term, 
with a coefficient that depends on Inm^. It will be seen in the next section that 
this is canceled by a counterterm insertion. There remain finite terms depending on 
the fermion mass that must be treated using dispersion relations when the internal 
fermions are quarks. 



6.3.2 Counterterm Insertions 

The (9(iVja^) corrections coming from 1-loop counterterm insertions into 1-loop dia- 
grams for the transverse part of the photon self-energy may be calculated to be 



It was checked that the longitudinal form factor vanishes which involves, once 
again, an intricate interplay between diagrams containing VF-boson 2-point countert- 
erms (0), scalar counterterms (|39| ) and vector-scalar mixing counterterms (^). 

By virtue of (BOF), the second term on the right-hand side of eq. ([75|) cancels the 



divergence in the fourth term of eq.(|7^) that depends on Inrrij. It may also be seen by 
substituting the expression for 5m f in eq.(0) into eq. (|76D that the remaining terms 
that depend on m/ in Il^Vy (0) are rendered finite by the counterterm insertions with 
the exception of the last. Moreover, in the on-shell renormalization scheme, there is 
an exact cancellation and these finite terms are eliminated completely. At 1-loop the 
bosonic and fermionic sectors of the theory are renormalized independently of one 
another. It is obviously a great convenience and simplification here to demand that 
the internal fermions are renormalized in the on-shell scheme. This will be done in 
the following but no such constraint will be imposed on the bosonic sector. 

Combining eq. ([75|) and eq.(|T6D we define a new quantity, 11^^^'' (0) = n^y^^''(0) + 
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n77'^Ho), given by 



2 \ 2 . 



x(!i^(.M|r-r(5-„)r(2-^)r(^-i 



+ 



167r2 / 3 



X7r2 ft^f{m))'i 2 

/ 

^^23^V 4(n3-12n^ + 41n-34) ^\^[^ 
+ VieTr^J n(n-3)(n-5) T 2^ T 2; 

xJ^QjK)"-^ (77) 
/ 

Of the two remaining terms in eq. (|77D that depend on m/, the first corresponds 
to weak corrections to photon-fermion vertex that have their origin in the diagram 
of Fig.^. The second comes from the pure QED diagrams, Fig.|^h&i, and their 
associated fermion mass counterterms. When expanded about n = 4 its leading 
logarithms reproduce the well-known result of Jost and Luttinger . Both sets can 
be treated via the dispersion relation trick, eg. . This requires that the diagrams 
be evaluated at some high q^. Fig-^e at high contains two unrelated scales and 
cannot be treated by the techniques used so far. In principle it can be obtained in 
closed analytic form from results given by Scharf and Tausk |2^ but it is neither 
compact nor illuminating and probably best obtained numerically. 

The pure QED diagrams of Fig.^&i are exactly calculable for Ig^l ^ m^. The 



24 



result is 



n;(2QED)^^2^ 



2 \ 2 

9 \ „A / ;-2\n-4: 



X <^ 8{n^ - 7n+ 16) 



+ 24- 



[n 



r(2-§)'r(§)V(§-2 

T{n)T{n - 1) 

'■-An + 8) r(4-r^)r(f)'r(5-2; 



n - l)(n - 4)' 



r(^-2 



(7^ 



Broadhurst et al. have given an expression for the subtracted photon vacuum 
polarization at general q^. The high-energy limit can be obtained by applying analytic 
continuation relations for the hypergeometric functions, 2F1 and 3F2 that appear in 
their result. 

Despite appearances, the expression on the right hand side of eq.([78D has only 
a simple pole with a constant coefficient at n = 4 that can be canceled by local 
counterterms. The leading logarithmic expressions can be found in ref. |30| , section 
8-4-4] where the authors invite the "foolhardy reader" to check that the finite parts 
are transverse. Here we have gone further and demonstrated this property in the 
exact result. 

It can be checked that the divergent part of ([78|) is identical to that of the last term 



in eq.(|77|) so that that the difference Rell^^ (g^) — Il^-y (0) is finite. This provides 
yet another useful check of various aspects of the calculation. 



6.4 Charge renormalization in a general scheme 

All the ingredients are now in place to write down the complete set of 0{NfQ^) cor- 
rections to the Thomson scattering matrix element and from it obtain an expression, 
extending eq.(^), for the physical electromagnetic charge of a fermion in terms of 
the renormalized parameters of the theory. 

Gathering all the parts together gives our main result 



^ \ ^2 ^^^ ^ ce Ml 



2 "ij I 2' 

So T Cfl 



g " 9' 



3 s ^.^^ + c^.^L^ (79) 

g 9 9 9 



where is given in eq.(^, ^z^i^) in eq.(0) and U^^^{0) in eq.(|3D. 
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7 0{Nfa^) wavefunction counterterms 



Up to this point we have been pursuing the expression for the physical matrix ele- 
ment for Thomson scattering and, to that end, terms that do not contribute to the 
final result have often been discarded. Although the final result does not depend 
on the 0{Nfa^) wavefunction renormalization counterterms, 5Z'^ and 5Z^^\ the 
Green's functions that have been calculated in the foregoing allow relations to be 
determined between them and dg^"^^ and Sg'^'^\ In all cases the leading divergences of 
these counterterms, i.e. those corresponding to a double pole at n = 4, are indepen- 
dent of renormalization scheme but subleading and finite parts will depend on which 
renormalization scheme has been chosen. 

When the 0{Nfa'^) diagrams contributing to Z-'-f mixing ( |66|) are combined with 
the counterterm insertions (|69D and the counterterms the result must be finite in 
any scheme. One thus obtains 



\ 9 9 9' 9' J 9 ^'2 2 

Turning to the 0{Nfa'^) corrections to the photon vertex it similarly follows that 
when the contributions from diagrams, (p5|), counterterm insertions, (p5D, and pure 
counterterms (^) are added together the result is finite. The part of the vertex 
proportional to yields the condition 

\n-4 , ^ 5g^^f^ ( g^ 



2 \ 2 




167r2y' 


n 






Ml 


Vl67r2 



r(4 - „)r (2 - |j r (^j + ^ 1 M4.)-TW 



^ V2 9 Vl67r2 



+ 3 • —■ — + — } o^^w + o^^B = finite 

\ 9 9 9 9 / 9 g' 2 2 

This expression differs from the previous one but is consistent with it because of 
the finiteness of the combinations ( p]| ) and (pT]). This consistency is a stringent check 
of very many aspects of the calculation. 
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The part of the the photon vertex proportional to Q leads to the condition 



. 6Z):' - 3 — .— = finite. (82) 



Obviously eq.s(0) and ( p2D can be used to obtain a condition for the combination 
of counterterms 

9 2'' c/ V2 9 J 9 9 

which is the OlNfO^) VT-fermion coupling counterterm and occurs, for example, in the 
calculation of corrections to the muon lifetime. Indeed this fact has been exploited 
as an extremely stringent cross check of both the results of this paper and of the 



calculation 0{Nfa ) corrections to the muon lifetime |^ 



The corrections to photon self-energy calculated in section |6.3| yield an indepen- 
dent condition on the the 0{Nfa'^) counterterms 

s^^zg) + 46 Z^^^ + n;(2) (0) = finite (83) 
where H^^ (0) is given in eg. (|77D . 



8 Summary 

The complete (9(A^/q;^) renormalization of the electromagnetic charge in the Standard 
Model using a general renormalization scheme has been presented. This represents 
the first practical calculation in which the full structure of the 2-loop renormalization 
has been confronted and lays the groundwork for future calculations of this type. The 
results have already been exploited in the calculation of the 0{Nfa^) corrections to 
the muon lifetime. 

The Z-7 mixing adds considerably to the overall complexity and number of Feyn- 
man diagrams that must be considered. Contributions coming from the insertion of 
1-loop counterterms in 1-loop diagrams constitute a significant portion of the calcu- 
lation due partly to the appearance at this order of counterterms that mix vector 
bosons with scalars. 

All integrals were performed exactly in dimensional regularization without ex- 
panding ine = 2 — n/2 yet in many cases they produced new and remarkably simple 
expressions that display the full analytic structure of the results. 

The role of wavefunction counterterms was investigated. It was found that the 
0{a) wavefunction counterterms, SZ^^\ had to be included in a manner consistent 
with Ward identities but that the 0{Nfa'^) counterterms, 6Z^'^\ could be neglected 
since they cancel in the final physical result. The price for this is that one must deal 
with divergent Greens functions in intermediate steps. 

A large number of internal consistency checks were performed in the course of the 
calculation in order to ensure the correctness of the results presented here. 
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